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Abstract 

Recent studies have pointed out the importance of 
the basal friction on the dynamics of granular flows. 
We present experimental results on the influence of 
the roughness of the inclined plane on the dynam- 
ics of a monodisperse dry granular flow. We found 
experimentally that it exists a maximum of the fric- 
tion for a given relative roughness. This maximum 
is shown to be independent of the angle of the slope. 
This behavior is observed for four planes with differ- 
ent bump sizes (given by the size of the beads glued 
on the plane) from 200 fim to 2 mm. The relative 
roughness corresponding to the maximum of the fric- 
tion can be predicted with a geometrical model of 
stability of one single bead on the plane. The main 
parameters are the size of the bumps and the size of 
the flowing beads. In order to obtain a higher preci- 
sion, the model also takes into account of the spacing 
between the bumps of the rough plane. Experimen- 
tal results and model are in good agreement for all 
the planes we studied. Other parameters, like the 
sphericity of the beads, or irregularities in the thick- 
ness of the layer of glued particles, are shown to be 
of influence on the friction. 



1 Introduction 

The flow of dense granular matter on inclined planes 
is often encountered in engineering applications in- 
volving the transport of materials such as minerals 



and cereals. It is also common in geophysical situa- 
tions where rock avalanches, landslides and pyroclas- 
tic and debris flows are natural events consisting in 
large-scale flows of grains Many chute flow 

experiments have been carried out and different con- 
figurations have been used, changing the boundary 
conditions from smooth E] to rough [§] , and using 
different kinds of materials |H1 El GUI without any 
systematic investigation. 

Literature on gravitational dense granular flows on 
inclined planes shows three regimes of motion de- 
pending on the angle of the slope. For small an- 
gles, the flow is decelerated; for intermediate an- 
gles, the flow reaches a constant mean velocity; and 
for large angles of inclination, the flow is constantly 
accelerated ^TJ . Recent studies have pointed out 
the importance of the basal friction on the dynam- 
ics of granular flow. They provide a model for 
slow frictional flows based on averaging equations (S* 
Venant equations') |12|. This model has been succes- 
fully tested on non-stationnary flows and is very 
promising for the dynamics of granular flow. Never- 
theless, one of the main parameter of the model is the 
basal friction. This friction is defined by the tangen- 
tial stress at the base divided by the normal stress at 
the base of the flow. In this paper, the word "friction" 
refers then to this global bulk effective friction of the 
flowing granular material down the slope and not to 
the local contact friction between two beads. This 
friction represents the interaction between the rough 
plane and the flowing layer. In this paper, we will 
not measure this friction but interpret experimental 
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results in term of increase or decrease of the friction. 
To our knowledge, no study has varied the roughness 
systematically in order to relate basal friction and ge- 
ometrical characteristics of the substratum. In fact, 
most of these studies were made with flowing beads of 
the same size ^3] or sometimes greater ^2] than the 
beads glued on the plane. It can be intuited that the 
most important parameter for understanding these 
flows will be the ratio between the size of the flow- 
ing beads and the size of the bumps (glued beads). 
Now, most industrial or natural flows do not prop- 
agate on a roughness composed of the same size of 
particles. Moreover, these flows are often polydis- 
perse in size. For example, geophysical events in- 
volve complex heterogeneous granular material made 
of particles of various size ranging from micrometers 
to meters, and often mixed with fluids. In such poly- 
dispersed material, it is well known that segregation 
occurs and that the large particles rise up to the free 
surface and go forth to the front ^JElEl- So, due 
to the phenomenon of segregation, the particles in 
contact with the substratum are different during the 
flow, involving different relative roughnesses. 

It also appears interesting to compare the motion 
of a granular mass down a rough inclined plane with 
the motion of a single bead down a rough plane. 
Many authors studied the trajectories of a single par- 
ticle moving down a rough surface. Unfortunately, 
most studies are limited to the case of a single par- 
ticle greater than the diameter of the glued beads. 
It has been found that the trajectories of the grains 
are very dependent on their size ratio, defined by the 
diameter of the flowing bead divided by the diameter 
of the glued beads. In this case too, it points out 
the fact that the relative roughness matters for the 
motion. A relative large grain has a greater probabil- 
ity to reach the bottom of an inclined surface than a 
relative small grain. Experimental 1171 IT51 119) . 
theoretical QH d EE E2 EH EH and numerical 
|rr7irrTfll!^n^ 1^1771 125] results show that the same 
three regimes of motion than those existing for granu- 
lar flows (decelerated flow, flow with a constant mean 
velocity, and accelerated flow) can be observed for 
one single bead. All these regimes depend mainly 
on the inclination of the plane and on the relative 
size of the bead. The geometry of the rough plane 



(random spacing 1201 EE]> regular spacing EE]) and 
the coefficient of restitution EH EE] play also a very 
important role on the motion of the particle. First, 
the angle for which there is a transition between the 
decelerated regime and the steady-state regime de- 
creases when the size ratio increases EE]- Second, in 
the decelerated regime, the characteristic length of 
trapping increases with the size ratio and with the 
angle of inclination |19j . Third, in the regime where 
the bead reaches a constant velocity, this velocity in- 
creases with the size ratio, and with the angle of in- 
clination EniEi), 28 . These results appear very inter- 
esting for modelling and understanding monodisperse 
flows down a rough inclined plane by considering one 
single bead on a plane. 

In this paper, we study several systems, constitued 
of different sizes of beads flowing on different planes, 
on which one layer of beads with a given diameter has 
been glued. It appears from these experiments that 
the relative roughness, defined as the ratio between 
the size of flowing beads and the size of glued beads is 
the main parameter. In order to study the influence 
of the relative roughness on a large range of value, 
two sets of experiments are performed: 

• varying the size of flowing beads on a given plane 

• varying the size of the beads glued onto the plane 

This study allows to determine for a large range of 
value (0.1 to 10) the influence of the relative rough- 
ness of a plane on the dynamics of a monodisperse 
dry granular flow. 

The three types of experiments and the measure- 
ment methods are presented in part 2. In part 3, 
we show experimentally the existence of a minimum 
of the spreading of a mass of granular matter, 
corresponding to a particular relative roughness. 
This minimum corresponds also to a minimum of 
the velocity of the flow and to a maximum of the 
thickness of the deposit. In part 4, we study the 
dependence of this extremum with few parameters. 
A model of stability predicts this extremum in 
part 5. The influence of some other experimental 
parameters is given in part 6. 
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2 Experimental configurations 

2.1 Experimental device 

The experimental device (fig. ^) consists in a 2 m 
long and 60 cm wide rough plane whose inclination 
(9) can be controlled precisely (less than 0.1° of er- 
ror). The particles are glass beads. We sieve the 
beads to obtain different ranges between f 50 fim and 
2mm, each range of beads is considered as monodis- 
perse, the width of the distribution being controlled 
by image analysis. These beads are used for flows and 
for rough planes (tab. [IJ. We avoid cohesive effect by 
using beads which have a size greater than f50 fim. 
A criterion of sphericity is taken as d max / d m in < 1.1, 
where d m ax and d m in are the two diameters of the el- 
liptic projection of the particle. The percentage of 
these non-spherical beads is less than 5 per cent in 
all batches of beads (tab. The plane is made 

rough by gluing one layer of particles onto the plane. 
For each plane (tab. we measure the compactness 
(c) which is equal to the ratio between the projected 
area occupied by all the beads and the total area. 
The compactness is an average on fifteen measure- 
ments each taken on a rectangle of a few hundred 
mm 2 area and a mean spacing is deduced. Error is 
given by the standard deviation of the fifteen mea- 
surements. We can also directly determine from the 
images of the plane, a value for the mean spacing 
which is the mean of all the spacings between the 
beads of the plane. With these results, we calculate 
the standard deviation obtained from the distribution 
of the spacing. The two values of the mean spacing 
are very close (tabic [5J • 

On each plane, three types of experiments are car- 
ried out (fig. HJ. The first set of experiments con- 
sists in the instantaneous release of a fixed volume of 
beads placed on the plane. The volume is placed in 
a hemisphere which is instantaneously removed. The 
released material flows down the slope, spreads and 
ultimately stops leaving a tear-shaped deposit, char- 
acterized by its length L and its width W. In order to 
control precisely the amount of the material present 
in the cap, the mass of beads poured into the cap is 
weighed before each run, and is equal to 176 g. In the 
second set of experiments, we study, for each range of 



beads and on each plane, the evolution of the thick- 
ness of the deposit (h stop ), left by a wide uniform 
steady state flow, for several inclinations of the slope 
(9). Third, we study, the velocity (u) of a stationary 
non uniform flow, produced by a constant flux input. 
Subsequently, we define the following parameters : 

• d : the diameter of flowing beads 

• A: the diameter of beads glued on the plane 

• c : the compactness of the layer of beads glued 
on the plane 

• e : the mean gap between two beads of the plane 
calculated from the compactness assuming a tri- 



angular lattice e = 




• 9 : the angle of inclination 

2.2 Measurement methods 

The thickness of the deposit is measured precisely by 
the deviation of the projection of a laser sheet. Two 
different pictures are compared: the first one is ob- 
tained with no beads on the plane (that is h s top = 0) , 
the second one with the deposit of a uniform steady 
flow. The deviation is proportional to its thickness 
{hstop)- The scale is determined by the deviation 
given by a plate whose thickness is known precisely. 
With this method, thickness is measured with less 
than 100 /im of error. 

For the measurements of the velocity of the non- 
uniform steady state flow, some marks are drawn 
along the plane, every ten centimeters. It allows us 
to measure the velocity on a video and also to verify 
if the flow is in a steady state. 

3 Existence of a maximum of 
the friction 

3.1 Results for the plane number 2 

In this part, all the experiments are carried out on the 
plane number 2. By keeping the same plane, and by 
varying the size of flowing beads, the flowing beads 
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Figure 1: Experimental set-up for granular flows on a rough plane: (a) Release of a volume contained in 
a cap; (b) Deposit left by a uniform steady state flow; (c) Free-side steady state flow. The plane is made 
rough by gluing beads on it. 



mean size {jJ/m) 


range 


(pm) 


mean size (/xm) 


range 


(pm) 


150 


140 - 


160 


580 


560 - 


600 


180 


160 - 


200 


615 


600 - 


630 


225 


200 - 


250 


655 


610 - 


700 


275 


250 - 


300 


670 


630 - 


710 


300 


290- 


310 


780 


710 - 


850 


327 


300 - 


355 


925 


850 - 


1000 


377 


355 - 


400 


1400 


1250 - 


1550 


425 


400 - 


450 


1850 


1700 - 


2000 


475 


450 - 


500 


2000 


1900 - 


2100 


530 


500 - 


560 


2150 


2000 - 


2300 








5000 


4900 - 


5100 



Table 1: Mean sizes of beads and corresponding ranges of sizes 



plane number 


1 


2 


3 


4 


diameter of glued beads (A) 


225 fim 


425 fim 


655 \im 


2 mm 


mean spacing (e) 


46 fim 


121 fim 


160 fim 


393 fim 


standard deviation of the mean spacing 


41 fim 


79 fim 


129 fim 


376 fim 


compactness (c) 


0.63±0.03 


0.55±0.04 


0.57±0.03 


0.65±0.05 


spacing deduced from c (e) 


44 fim 


120 fim 


171 fim 


362 fim 



Table 2: Rough planes 
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are submitted to a variable relative roughness. In 
fact, one plane can be rough for some beads, and 
smooth for others. The couple (plane, flowing beads) 
allows us to define the relative roughness of our sys- 
tem. Another mean to vary this relative roughness 
would be to vary the size of glued beads and not the 
size of flowing beads, but it is less precise because it 
is difficult to control the spacing between the glued 
beads with accuracy. 

First, the experiments of the release of a constant 
volume are done to study the length of spreading 
(L) versus the diameter of beads flowing (d). These 
experiments arc carried out for different inclinations 
(0). Results are presented on figureEJ For an angle of 
inclination greater than 28° , the measured lengths arc 
underestimated for beads which have a size greater 
than 327 pm (425 fim, 475/im, 530 fim), because 
some beads do not stay on the plane. For these sizes, 
the percentage of loosing beads is respectively equal 
to 3, 8, 13 %. For an angle of inclination smaller 
than 28°, all the beads stay on the plane. The length 
presents a minimum for a given size of beads, named 
d c (fig. 0). This diameter is independent of the an- 
gle of inclination. But, the minimum of the length 
of the spreading is more enhanced for large angles of 
inclination. 

On the previous experiments, we also measure the 
maximal width (W) of the deposit, for different sizes 
of beads (fig. [3J|. The maximal width is independent 
of the size of the beads flowing. 

Second, experiments of wide flows give the ability 
to study the thickness of the deposit (h s t op ) left by a 
uniform steady state flow. For these experiments, the 
inclination and the size of flowing beads vary. Figure 
0] presents the curves of h s t op as a function of 0. The 
curves are classified, all over the range of inclinations 
we study, and h stop is greater for d — d c than for 
other sizes. Figure presents h st0 p as a function of d 
for different angles of inclination. For d — d ci h stop is 
greater than for other sizes, and that for every angle 
of inclination. 

Third, velocity (u) is measured on steady state 
free-side flows (fig. EJ). As in the previous experi- 
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Figure 2: Length of the deposit as a function of d 
for different angles of inclination: (■) =28,3°; (k) 
61=25,7°; (•) 0=22,8°; (♦) 0=18,4°. There is a diam- 
eter d c , for which the length presents a minimum. d c 
is independent of 9 (plane 2). 
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Figure 3: The maximal width of the deposit is inde- 
pendent of d (0=25,7°, plane 2). 
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ments, the behavior of the flow presents a singularity 
for d = d c : the velocity is minimum for this value. 



24 26 28 30 32 34 36 
angle of inclination 6 (degrees) 

Figure 4: h stop for different sizes of beads on plane 
2: (•) d=150 fim; (a) d=225 [im\ (♦) d=275 ^m; 
(A) d=327 ^m; (□) d=425 fim; (o) d=530 /xm. For 
d = d c , the /istop curve is maximum (c? c = 275^to). 
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Figure 6: Velocity as a function of d. For d = d c , it 
exists a minimum of the velocity (0=25.7°, plane 2). 
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Figure 5: h st0 p for different angles of inclination: (•) 
6>=27°; (■) 0=28°; (A) 6»=28.3°; (♦) 6>=30°. For 
d = d c , h s top is maximum (plane 2). 



We notice that the curves L, h stop and u as a function 
of the diameter of the flowing beads d are not sym- 
metrical in relation to d = d c . In these three types 
of experiments, we point out that it exists a value of 
the diameter of flowing beads, named d c for which : 

• the thickness of the deposit (h st0 p) is maximum 

• the length of the deposit (L) is minimum 

• the velocity of the non-uniform steady state flow 
(u) is minimum 

The experiments show that the diameter d c is inde- 
pendent of the angle of inclination of the slope 8. 
These results can be linked to the basal friction of 
the flow. For d < d c , the friction increases with the 
diameter of beads flowing. For d > d c , the friction de- 
creases when the diameter of beads flowing increases, 
and there is a maximum of this friction for d = d c . 

To found the same extremum on these 3 parame- 
ters could be predictable, because the three parame- 
ters, L, h s top and u are linked. We approximate the 
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angle of inclination 9 (degrees) 

Figure 7: h s t op for different sizes of large flowing 
beads: (□) d=275 (j.m; (o) d=530 fim; (♦) d=l mm; 
(•) d=2 mm; (a) d=5 mm (plane 2). The curvature 
of h s top tends to and 9\ decreases for large beads. 

surface of the deposit by a rectangle and assume that 
the volume (V) is equal to: 

V = WLh stop (1) 

and as the volume of flowing beads is constant, de- 
duce that L and h stop vary as the inverse of each 
other because W is constant (fig. [3}. On the other 
hand, if we suppose that results given by Pouliqucn 
|12j are verified in our configuration, the velocity of 
a steady state flow is a function of h s top and of the 
height of the flow (h): 



\fgh h stop 

with (3 a constant. A maximum for h stop at d = d c 
is then equivalent to a minimum of the velocity u at 
the same value of d. 

We could interpret the existence of the minimum 
of the spreading for the two borderline cases: d^S> X 
and d -C A. 

For d 3> A, the length L and the velocity u in- 
crease with d and with the angle of inclination 9 (fig 



13 EJ In fact, beads are greater than A and than 
the spacing between the glued beads of the plane (e), 
the relative roughness is small and beads flow eas- 
ily. Flows are submitted to decreasing friction when 
d increases. This behavior is in accordance with the 
results of previous studies on the motion of a large 
single bead whose velocity ^| [221 EH or character- 
istic length of trapping 19 increases with the size 
ratio and with the angle of inclination 9. The range 
of angles of inclination [0i;#2] (for 9 — 9\, h s t op —* 
oo, for 9 = 92, hstop— 0) for which a steady state flow 
happens, decreases when d increases . For greatest d 
(figEJl, this range of angles of inclination is very small 
(d=l mm) and tends to (d=2 mm, d=5 mm). At 
the same time, both 9\ and O2 decrease. There is 
a limit for which we can only define one angle. For 
an angle greater than this angle, no deposit stays on 
the plane. For an angle of inclination lower than this 
angle, no steady state flow is observed. When d in- 
creases, the value of this limit angle decreases, beads 
flow easily. The global shift of h st op curve toward 
small angles when d increases could be related to the 
decrease of the angle separating steady motion from 
the decelerating motion of one single bead on a rough 
plane [2B]> although, the values of this angle are very 
small (approximatively 9° for a size ratio equal to 1) 
compared to 9\ and #2. 

For d -C A, the length L increases when d de- 
creases, which is more surprising. An explanation 
could be that the holes filled with small beads reduce 
the roughness. 

But this interpretation of the two borderline cases 
do not allow us to explain all our results because the 
evolution of the length with the diameter seems con- 
tinuous and does not only show a transition below a 
critical size which fills holes. Consequently the varia- 
tions of the friction with d is not due only to the filling 
or not of the holes of the plane. We will discuss more 
precisely this point in part 4. 

3.2 Results for different planes 

In this part, we study the behavior of beads on the 
different planes, presented in part 1. We measure 
for each plane and for different ranges of beads, the 
length L of the spreading of a constant mass and 
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hstop- For all these planes, there is a diameter of 
beads which presents a minimum of length L and a 
maximum of h s t op . We compare the diameter of these 
beads (d c ) with the diameter of beads glued on our 
plane (A) (fig. EJl. 

The beads which correspond to the minimum of 
length are not always the same. This minimum is 
then not an experimental artefact due to some char- 
acteristics of a batch of beads. 

d c is approximatively equal to A/2 with little 
deviations from this value. But, in all of our 
experiments, the variation of d is discontinuous, and 
we can not determine with accuracy the value of d c 
for which the friction is really maximum. 




diameter of glued beads A, (|J.m) 

Figure 8: d c obtained on different planes are close 
to A/2 (dotted line): (•) experimental results; (□) 
model with the experimental parameters values. 
Curves represent predictions of the model taking into 
account of different compactnesses (c). The model 
predicts with accuracy the value of d c . For the model, 
error bars (O)come from the fact that c and A are 
not well-determined. For the experiments, d c is not 
known precisely because of the discontinuity of the 
values of d (-). The error bars (axis A) for the model 
and the experiments are the same, taking into ac- 
count of the range of sizes of the glued beads. 



4 Implications on the rough- 
ness of filling the holes 

It appears also interesting to understand the impact 
on the friction of filling the holes of the plane by small 
beads. We used the plane number 4 made by gluing 
2mm beads, so it is possible to study a wider range of 
small d without using cohesive small beads less than 
150 fj,m. Small beads, for which the ratio d/X is very 
small, fill the holes of the plane. For example, for 225 
pjm beads, one hole of the plane contains approxima- 
tively 30 beads. The measurements of h stop and L 
are carried out on plane number 4 and only some re- 
sults of h s to P are presented on figure El The behavior 
is similar to the one observed previously, with h s t op 
maximum and L minimum for d c =925 /im. 

If we compare the thickness of the deposit in term 
of number of beads (h st0 p/d), we observe that (fig. 
EH: 

• d c could be a local extremum (the friction could 
be greater for the smallest beads) 

• the curves corresponding to <i=225 fim and 
d=327 fim are the same. 

For the smallest beads, the value of h s t op /d seems 
to tend to an asymptotic value. 

This behavior was not observed for the other 
planes, because the ranges of the flowing beads com- 
pare to the glued one, were not wide enough. 

For the smaller beads, all the holes of the plane are 
filled, and the basal plane is a mixture of large beads 
(2mm) and small beads (225 \im or 327 fj,m). One 
possibility to interpret the fact that h st0 p/dior d=225 
\xm and for d=327 fim are equal is to consider that 
the friction comes principally from the filled holes and 
that the contribution of the summit of large glued 
beads to the friction is negligible. We can suppose 
that the beads flow on a plane which has a rough- 
ness given by the size of beads organized with a 3D 
random compactness in the holes. In that case, the 
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Figure 9: h stop on the plane number 4 (2mm glued 
beads) for different angles of inclination: (•) 9=24°; 
(■) 6>=26°; (♦) 6=27°; (▲) (9=28°; (+) 6»=29°; (▼) 
0=30°. There is a diameter d c =925 \im for which 
^stop is maximum. 



Figure 10: h s top/d on the plane number 4 (2mm glued 
beads) for different angles of inclination: (•) 0=24°; 
(■) 0=26 O ; (♦) 9=27°; (A) (9=28°; (+) (9=29°; (T) 
9=30°. There is a diameter <i c =925 fim for which 
hstop/ d is locally maximum. 



behavior is independent of the size of the beads ex- 
plaining the similarity between the curves for 225 fim 
and 327 [im. 

Because of the filling of the holes, we also suppose 
that 225 fim flowing on 2 mm rough plane could be 
similar to 225 /im flowing on 225 /im rough plane. It 
is so interesting to compare 225 /im beads flowing on 

• a 225 fim plane (c = 0.63), plane 1 

• a 2mm plane (c = 0.65), plane 4 

In the second case, the beads flow on a plane con- 
stituted of 225 /im not glued small beads and very 
smooth rounded surfaces (top of 2 mm beads). In- 
side the holes, the lattice of 225 fim beads is randomly 
distributed with a compactness between 0.5 and 0.7. 
This compactness is also the compactness of a plane 
section of a random three dimensional piling. 

The h s top curves are very close each other except at 
small angles (fig. II lfl . From these results, we assume 
that 225 \xm beads flowing on a 2 mm plane could 



be equivalent to 225 \xm beads flowing on a 225 fim 
plane which have a compactness equal to the com- 
pactness of a random lattice. This result is surprising 
because it means that there is almost no influence of 
the smooth part of large beads tops, at least when 
the deposit is thin. For small angles of inclination, 
the difference could be explain by the difference in 
the values of the compactness. 

Notice that for 225 /zm beads flowing on a 2mm 
plane, filling the holes corresponds to introduce 
rough surface between the large beads. But for 
2mm flowing beads, filling the holes with 225 fim 
corresponds to a smoothing of the plane. 
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Figure 11: h stop as a function of 9 for the same 
225 fim beads flowing, on different planes: (•) plane 
number 1 (A = 225^m); (□) plane number 4 (A = 
2mm). The curves represent the data of h s to P fitted 
by the function (J3J) and the two associated parameters 
9\ and 8 2 - 



5 Model of stability for the ex- 
istence of the maximum 

5.1 Stability of a single bead on a 
rough plane 

We have seen in the previous part, that the param- 
eters, h s top, L and u, all present an extremum for 
the same value of d c . We consider, in this model, a 
rough plane defined by two parameters: the diame- 
ter of beads glued A, and the mean spacing between 
two beads glued on the plane e. The model con- 
sists in determining the angle corresponding to the 
limit of stability of one single bead put on the plane. 
The criterion of choice of the exact location of this 
bead is such that the fall of this particular bead put 
h s top equal to 0. Among all the possible locations, we 
choose the more stable one, which would correspond 
to the last bead to fall before h stop = in the case of 
a whole layer of beads covering the plane. According 



to this criterion, we calculate an angle corresponding 
to the limit where there would be no more deposit in 
the configuration of our experiments. 

This angle can be compared to another angle, 9 2 
obtained from the h stop versus 9 experimental curves. 
6>2 is defined by the angle for which h stop (9)=0 and 
would then correspond to our criterion previously ex- 
plained. In fact, from each curve h stop as a function 
of 9, we can deduce the two parameters, 9\ and 9 2 
from the curves (fig. Illfl . For 9 < 6\, no uniform 
steady flow can be observed (h stop — oo). For 9 > 9 2 
no deposit stays on the plane (h stop = 0) ^2]- We 
can also determine more precisely the values of 9 2 
by using the whole set of experimental data of h s to P 
fitted by the following function [T2| : 



tan9 — tan9i + (tan9 2 — tan9\) exp 



*'StO P 

~Ld~ 



(3) 



We also experimentally determine that the level 
h s to P = corresponds to the plane tangential to the 
top of the beads glued on the plane, because there is 
no deviation of the laser projection when the holes 
are filled with very small beads. That means that a 
plane where only the holes are filled with small beads 
will be considered as having no deposit on it with our 
measurement method. From these observations, we 
can define a criterion for the choice of the bead of 
which we will consider the stability: the bead cho- 
sen is the lowest bead which top is above the level 
h stop =0. Consequently, if the bead is small enough 
to fit between the glued beads of the plane, the lowest 
bead is the bead, which top is tangential to the level 

h stop =0. (fig. [13 

In our model, we determine the minimum angle 
for the chosen bead to fall down. This angle is 
named the angle of stability. Our model is then 
based on a phenomenon of start because we consider 
the stability of a bead, which is initially at rest. 
The angle of stability, we determine, is the angle 
for which this bead becomes unstable. However, 
9 2 is experimentally the angle for which no deposit 
stays on the plane after a flow. Consequently, 
the measurements of h stop are determined by a 
phenomenon of stop. But neglecting the effect of the 
inerty in the process of stop, the angle of stability is 
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then equal to 82- 

We first explain our model of stability in a two- 
dimensional case, the same arguments being valid 
in the three-dimensional case. The results presented 
in the figures iJSJ 1141) are from the three-dimensional 
case calculations. 



5.1.1 The two-dimensional case 

The bead which stability we consider is tangential 
to one glued bead of the rough plane, and its top is 
above or tangential to the level h stop =0. We define 
the diameter (dh), which is the diameter of the bead 
tangential to the two beads of the rough plane, and 
tangential to the level h st o P —0 (fig- El case b). For 
beads greater than dh, the chosen bead is above the 
level h s to P =0, and tangential to the two beads of the 
plane (fig. El case a )- For beads smaller than dh, the 
chosen bead, is tangential to the level h s to P =0 and to 
one bead of the rough plane, the holes are filled with 
small beads (fig. El case c )- 

Once the location of the bead chosen, the criterion 
of stability is very simple. The bead is submitted to 
2 forces, its weight, named P, and the contact force 
between the bead and one bead of the plane, named 
R (fig. I12|) . The direction of R depends on the angle 
of inclination of the plane 8. We consider the angle 
between these two forces, 4>. If this angle <fi is equal 
to 0, the bead is unstable, else the bead is stable. 

For the cases presented in figure we calculate 
the angle of inclination for which the bead, initially at 
rest, is unbalanced. This angle is equal to 82 because 
for an angle lower than the angle of stability of this 
bead, h st0 p is not equal to 0, and for an angle of 
inclination greater than this angle of stability, the 
bead falls leading to h stop =0. 

With this choice of the unstable bead, and with 
this criterion of stability, it appears that the angle 
of stability is maximum for the bead whose centre is 
the lowest; i.e. for the bead with the diameter dh- 
We realize with this schema that the maximum of 
the friction comes from a pure geometrical reason. 




(a) (b) (c) 



Figure 12: Criterion of stability in different cases: (a) 
d > dh', (b) d = dh and (c) d < dh- The dotted line 
represents the level h stop = 0. The stability of the 
bead is maximum for d = dh- 

5.1.2 The three-dimensional case 

In that section, we consider the stability of the bead 
which is tangential to 2 beads of the rough plane, 
and which top is still above or tangential to the level 
h s top— 0. We assume that the plane is ordered with 
a triangular lattice (fig. I13fl . We calculate the di- 




(1) 

Figure 13: Representation of the roughness plane tri- 
angular lattice, with a spacing e. d m is the diameter 
maximum of the bead that can be placed in between 
the three beads of the rough plane. The direction (1) 
is for the bead which pass over the gap between two 
beads of the rough plane, the direction (2) is for the 
bead which get over one glued bead, dh is the diam- 
eter of the bead tangential to the three glued beads 
and to the level h stop =0. 

ameter maximum (d m ) of a bead, lodging between 
the three beads glued onto the plane. We also define 
the diameter (dh) of the bead which is tangential to 
the three beads of the rough plane and to the level 
h s top = 0. This diameter is equivalent to dh defined 
in the two-dimensional case, but their values are not 
the same. 

The criterion of stability is the same, but we could 
determine a whole range of angles of stability. The 
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smallest, is for the bead which passes over the gap d very close to d m , with d, 
between two beads of the rough plane (fig. El direc- 
tion 1), the highest is for the bead which gets over 
one bead of the rough plane (fig. El direction 2). In 
the following, we only consider the smallest angle of 
stability (fig. El direction 1). 



_ 2(A + £ ) 



(4) 



The choice of the bead is separated in several cases: 
for d > dh, the case is well-defined: the chosen bead 
is tangential to the three beads of the rough plane. 
The calculations, for the case d > A, have been made 
by Riguidel ^2]. Here, we obtain the same results 
for d > dh- For d < dh, we study two cases. For the 
largest d < dh, the chosen bead is tangential to two 
beads of the rough plane and to the level h stO p—0. For 
the smallest d, the bead is too small to be tangential 
to two beads and to the level h stO p=0- In that case, 
the beads, which are also smaller than d m , fill the 
holes between the large glued beads. Because there 
is a filling of the holes, we are not sure that the more 
unstable bead is one touching a large glued bead, or 
one being on the surface of the filled hole. At the 
limit d — >0, we expect that there is no more effect of 
the top of the large glued beads on the stability of 
the chosen small beads. We choose to consider the 
stability of a small bead which stays on the random 
loose packing of the others small beads. For that 
reason, the angle of stability would not vary with the 
size of the small beads for d — ->0. For the calculations 
involving the smallest beads, we thus consider the 
stability of a small bead tangential to three small 
beads and to the level h stO p=0 (being at the top of the 
random piling which fills the holes); the three small 
beads spacing being determined by the assumption 
of a random packing in the holes between the large 
glued beads. In all our calculations, we choose the 
compactness of the surface of a random packing equal 
to 0.57. We discuss this value of the compactness 
in part 5.2. The diameter for which, we pass from 
one calculation to another is not well defined in our 
model. Consequently for d < d h , the angle of stability 
will be taken equal to the maximum of the 2 angles 
given by the two cases. 

The curve is decreasing for large beads and increas- 
ing for small beads (fig. 1141) . For small beads, the 
intersection of the two calculations of 02 happens for 



Consequently the stability of beads smaller than d m 
is the same wathever their size are. 

The main result is that the angle of stability 02 
presents a maximum for dh' 



_ (A + e) 2 tt 
dh 3A ~ ^3 



(5) 



One more time, we notice that the maximum is due 
to a pure geometrical reason. The model predicts 
that dh is proportional to A but also that it depends 
on the spacing of the glued beads e. 

We have also calculated the angle of stability, 
considering the possible patterns fitting in the holes 
according to the size ratio between flowing and glued 
beads. But it appears that the results are not really 
different except that there are some discontinuities 
when passing from a pattern to another. The global 
behavior being the same, we choose not to represent 
the results. 



5.2 Comparison between the exper- 
imental results and results ob- 
tained from the model 

In figure ITT1 we compare, for one plane, experimen- 
tal results, and results obtained from the model. We 
have also compared the results for the others planes; 
the same following remarks being valid. The global 
behavior of experimental results and model are the 
same: for d < d m 02 is constant, for d m < d < dh, 
the angle 02 increases with d, and for d > dh, 02 de- 
creases. The curvature for experimental results and 
results predicted by the model seems to be the same. 
Taking into account of the discontinuity in the ex- 
perimental variation of d, the model predicts with 
accuracy the value of the extremum, and in the fol- 
lowing text, d c is put equal to dh- In contrast, for 
small diameter of the flowing beads, the value of 02 
given by the model is underestimate. In fact, we have 
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seen that the criterion of stability is very difficult to 
apply precisely for the small beads. The compactness 
of the top of the random piling we have choosen is 
equal to 0.57. However, all we know is that the value 
of the compactness is between 0.3 and 0.9 (see the 
end of this section), which respectively correspond 
to an angle of stability 82 equal to 90 j and 19,7 j. 
This range of values of the compactness is not pre- 
cise but we do not consider valuable to calculate ac- 
curately the fraction because the top of the random 
piling is probably not equivalent to a flat pattern of 
small glued beads with a certain spacing. Neverthe- 
less, the model shows that it exists a maximum and 
predicts with accuracy the value of d corresponding 
to the maximum (fig. |HJ. 

In figure 00 dh obtained from the model for all 
planes, are plotted as a function of the diameter of 
glued beads (A). The model predicts that dh is pro- 
portionnal to A, when no variations of the compact- 
ness. But the compactness changes from one plane 
to another explaining that the values of d c are not 
aligned. We expected that all the experimental points 
will be between the two lines corresponding to com- 
pactness 0.5 and 0.7, in accordance with the mea- 
sured compactnesses of the plane (tab. 2) which is 
approximatively verified. 

Errors bars for results from the model (fig. |SJ 
take into account of the uncertainty standard de- 
viations on the values of the compactness because 
the measurements could be not representative, the 
surface analysis being only equal to few thousand 
mm 2 . Moreover, we suppose in the model that only 
the mean value of the compactness matters and not 
the variability of the spacing around its mean value. 
Errors bars for the model take into account of the 
standard deviation of the direct measurements of the 
spacing, and suggest that this variability is impor- 
tant. Nevertheless, assuming values for A and e, the 
model predicts to our satisfaction the value of d c . 
However, molecular dynamic simulations on the the 
motion of a single particle on a bumpy inclined plane 
[281 |2*U| show that intoducing disorder has a similar 
influence on the motion of the ball as the introduction 
of a regular spacing. The velocity of the ball is not 
affected by the introduction of disorder to the line, 
the velocity can be approximated by the mean veloc- 



ity on a line with equally spaced balls: this spacing 
corresponds to the mean value of the disordered case. 

Because experimental reasons, all the planes, we 
made, have also a compactness between 2 values: 
Cmin ancl Cmax- Cmax is the compactness for which 
all the beads glued are in contact. This compactness 
is equal to c ma x = 7r/2\/3 = 0.9 for monodisperse 
size beads. The minimum of the compactness is 
determined by the manufacturing process. To make 
the planes, we put a large amount of beads on a 
sticky paper, we reduce the spacing between the 
beads by pressing on the beads and remove the 
beads not glued. No hole larger than a bead size 
stays empty. The minimum of the compactness is 
obtained when the spacing between three beads 
(d m ) is equal to A, leading to a compactness equal 
to c m i„=0.3. With these extreme compactnesses, 
the model predicts that the value of the diameter d c 
is between these two lines, and consequently, d c is 
always less than A (fig. |HJ). 

70 




20 

100 200 300 400 500 600 

diameter of the flowing beads d (urn) 



Figure 14: 82 as a function of d. (A) experimental 
results obtained on the plane number 2 (A=425 /im, 
c=0.55) (□) experimental results obtained by the fit 
(|3*|l. The curve represents the results given by the 
model with A=425 \xm and c = 0.55 
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6 Influence of others parame- 
ters 

In the previous part, we worked with rough planes 
made of one layer of spherical beads. It appears also 
interesting to study the influence of the quality of 
beads and planes. First, we will study the influence 
of the sphericity of glued and flowing beads, "non- 
spherical" beads batch is composed to approxima- 
tively of 60% of spherical beads, of 15% of splinters 
and of 25% of non-spherical beads, with the criterion 
of sphericity defined in part 2. Second, we will study 
the influence of a mutli-layered rough plane whose 
amount of layer is spatially variable. In figure 1151 
we study four different configurations for the length 
of the spreading (L) as a function of the diameter of 
the flowing beads d: 

• (1) spherical beads flowing on the plane number 
1 on which one layer of spherical beads are glued, 
named "good" plane 

• (2) non-spherical beads flowing on the "good" 
plane 

• (3) spherical beads flowing on a plane on which 
there are, depending on the location, approxima- 
tively one to three layers of 225 fim non-spherical 
beads, named "bad" plane 

• (4) non-spherical beads flowing on the "bad" 
plane. 

6.1 Influence of sphericity of beads 

In this part, we compare the experimental results ob- 
tained with the configurations 1 and 2 (or 3 and 4) . In 
each case, spherical and non spherical beads flow on 
the same plane. We notice that it exists, in these four 
cases, a diameter for which the length of spreading is 
minimum. The lengths obtained on the same plane 
are very similar although the lengths for non spheri- 
cal beads are 1.6 smaller than for spherical beads in 
the cases 3 and 4. The sphericity of the flowing beads 
seems to have just a little influence on the spreading. 
In fact, during the flow, because of the segregation, 
the splinters, which represents 15% of the volume, 



segregate at the surface and at the front, and they 
consequently are probably of little influence on the 
flow basal friction. 



6.2 Influence of a multilayered rough 
plane 

In this part, we compare the experimental results ob- 
tained with the configurations 1 and 3. Spherical 
beads flow on a "good" plane and on a "bad" plane. 
We see that the quality of the plane is of great impor- 
tance. If we compare the configurations 2 and 4, done 
with "non-spherical" beads, we see the same effect. 

First, the length of spreading is greater in the case 
of a "good" plane than in the case of a "bad" plane. 
This result may be explain by the fact that the pres- 
ence of splinters increases the roughness of the "bad" 
plane. Moreover, the steps between the layers cre- 
ate a barrier difficult to get over. That is why the 
lengths of flowing are really smaller in the case of a 
"bad" plane. 

Second, we observe, that in the case of a "good" 
plane, for a size of beads greater than a limit size 
(275 fim for these experiments), no deposit stays on 
the plane (length is arbitrary marked as 200 cm in 
that case). But in the case of a "bad" plane, there is 
a deposit for larger flowing sizes. The behavior of the 
beads, in the case of they are greater than d c is totally 
different. This behavior can be also interpreted as the 
"bad" plane being more rough than the "good" one. 

We notice that it exists always a diameter d c for 
which the length is minimum. Unfortunately, our 
data are not precise enough to conclude if the value 
d c is the same in the two cases. 

In conclusion, the sphericity of flowing beads has 
less influence, the global behavior is not affected by 
the quality of the bead. But the quality of the plane 
is of great importance, the main effect is that the 
relative roughness of a plane can be totally smooth 
(cases 1 and 2) or still rough (cases 3 and 4) for large 
flowing beads. 
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Figure 15: Length as a function of d: (♦) spherical 
beads flowing on a one-layer spherical beads plane 
(case 1); (•) non-spherical beads flowing on a one- 
layer spherical beads plane (case 2); (a) spherical 
beads flowing on a multi-layer plane (case 3); (■) 
non-spherical beads flowing on a multi-layer plane 
(case 4). L = 200cm represents experiments where 
no deposit stays on the plane (A=225 \im 9=25.7°). 

7 Conclusion 

This paper presents experiments of flows of various 
sizes of glass beads down rough inclined planes. Vary- 
ing the size appears equivalent to a variation of the 
relative roughness of the plane. We show the exis- 
tence of a size of beads (d c ) for which the length of 
the spreading and the velocity are minimum and the 
thickness of the deposit is maximum; these facts be- 
ing interpreted as a maximum of the friction. This 
behavior is robust, being observed for several planes 
of different roughnesses. We also study the influence 
of different parameters on the value of the diameter 
d c corresponding to the maximum of relative rough- 
ness. First, it is independent of the angle of inclina- 
tion. Second, it is mainly determined by the size A of 
the particles glued on the plane, but the compactness 
also has an important influence on it. When the com- 
pactness of the rough plane increases, the value of d c 



decreases. A simple geometrical model of stability of 
one bead allows to estimate the value of this diameter 
d c , by calculating the value of the angle of inclination 
above which no deposit stays on the plane. 

Most flows are more complex than the simple 
avalanche of monodispersed glass beads. The interac- 
tions between the plane and the flow will be difficult 
to understand because each type of particle experi- 
ences a different relative roughness. For example, in 
the case of a bidisperse granular flow, the existence 
of a minimum for d c leads probably to complex be- 
haviors depending on the positions of each sizes com- 
pared to d c combined to the segregation effect. Our 
results provide a good frame for the interpretation of 
the dynamics of the flow of a granular polydisperse 
matter down a rough plane. 
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